Ergodicity of Random— Matrix Theories: The Unitary Case 
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We prove ergodicity of unitary random-matrix theories 
by showing that the autocorrelation function with respect 
to energy or magnetic field strength of any observable van- 
ishes asymptotically. We do so using Efetov's supersymmetry 
method, a polar decomposition of the saddle-point manifold, 
and an asymptotic evaluation of the boundary terms gener- 
ated in this fashion. 
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Predictions of random-matrix theories (RMT) are ob- 
tained by averaging over the ensemble of random matri- 
ces. To test such a prediction, one compares it with data 
obtained from a single system (rather than from an en- 
semble of systems governed by different Hamiltonians). 
For the data set, the ensemble average is replaced by the 
running average (i.e., the average over energy E, mag- 
netic field strength B, or another variable). If the two 
averages are, in a physically plausible limit, almost al- 
ways equal irrespective of the observable at hand, RMT 
is said to be ergodic. 

It is generally believed that RMT is ergodic. How- 
ever, proofs have been given only for spectral fluctua- 
tions and some other special observables, see the review 
[]. In this Letter, we report on the first general proof 
of ergodicity for the Gaussian Unitary Ensemble (GUE), 
the ensemble of Hcrmitean matrices H^ v of dimension 
N where N — > oo. The matrix elements are Gaussian 
random variables with zero mean and a second moment 
given by H^H^^ = (A 2 / 'N)S \ IV * S M > v . Our proof is not 
restricted to this ensemble, however, and applies equally 
to all other non-Gaussian unitary ensembles specified in 
the paragraph below Eq. (^). This is because all these 
ensembles are known [] to possess identical local fluctu- 
ation properties. It is these fluctuation properties which 
are used in the proof. In the framework of the GUE, 
we consider the two most important cases: The running 
average extends either over energy or over magnetic field 
strength. 

The proof is based on the use of Efetov's supersymmet- 
ric generating functional [] and on a polar decomposition 
of the saddle-point manifold. The latter yields a series 
of boundary or Efetov-Wegner terms which we generate 
using Refs. []. The restriction to the unitary case is for 
technical reasons only. We feel certain that the method 
of proof works equally for the technically more demand- 
ing orthogonal and symplectic ensembles. Lack of space 
forces us to give the essential steps only. Details may be 
found in Ref. [1. 



In the standard approach to the problem [] one consid- 
ers an arbitrary observable F. To be specific, we assume 
that F is a function of energy E, so that F = F(E,H) 
where H denotes the Hamiltonian. We return to the case 
of the magnetic field strength at the end of the paper. We 
denote the ensemble average (the running average) by 
an overbar (an angular bracket, respectively). Explicitly, 
FjE) = J dfi(H)F(E, H) where dfi(H) denotes the GUE 
measure for integration over H (including the Gaussian 

weight factor), and (F(E)) = A" 1 dE'F{E',H) 
where A denotes the averaging interval, and where H 
is an arbitrarily chosen member of the random-matrix 
ensemble. To define ergodicity, one asks whether 



(F{E) 



(F(E)) 



for all F(E, H) 



(1) 



The arrow denotes a limit defined below. If this condition 
is met, the difference between the two averages tends 
almost always to zero, and the ensemble is said to be 
ergodic. 

The GUE spectrum is confined to the interval — 2 A < 
E < 2A, so that the average level spacing d ~ X/N. 
In the ergodicity test (P, one considers the limit where 
the running average (RA) is performed after the limit 
N — > oo has been taken. The RA extends over K lev- 
els where K — * oo and, thus, over an energy interval 
of length A = Kd which, for any finite K, is infinitesi- 
mally small in comparison with the range 4A of the spec- 
trum. The ensemble average F(E) depends on E only 
via the mean level density p(E) and is, on the scale Kd, 
therefore constant (local stationarity) . It follows that 
(F(E)} — F(E). This fact allows us to cast the ergodic- 

The 



ity condition in the form {F(E)) 2 - [F(E) 



0. 



term (F(E)) 2 can be written as a double RA taken over 
the function F(E 1 )F(E 2 ) = F(E U H)F(E 2 ,H). This 
function is symmetric in E\,E2, and, for E\,E 2 C A, 
depends only on \E\ — E%\. With the angular brackets 
denoting now this double RA, the ergodicity condition 
reduces to 



F(E 1 )F(E 2 ) - 







(2) 



Application of Slutsky's theorem [] shows that ergod- 
icity holds provided the autocorrelation function Cg — 
F(E 1 )F(E 2 ) - F{E X ) F(E 2 ) vanishes for large \E X - 
E 2 \/d. This is what we are going to show. 

In general, F(E, H) may itself be a fc-fold correlation 
function. Then, F(E, H) will primarily not depend on 
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E but on the k energy arguments ex, . . . , where A: is a 
positive integer. To implement the ergodicity limit, we 
write Cj = E + u>j with J^j w j = an d E = (1/fc) J2j e j- 
The ujj 's are defined on the scale d and are held fixed, and 
we consider F as a function of E only. Using standard 
notation and procedure [], we can cast F(Ei)F(E2) in 
the form 

F{E l )F{E 2 ) = J d A1 (T)exp(z7r[e/d]trg(Qr 3 ))5(Q) . 

(3) 



to vanish asymptotically for large e/d, except for a term 
which cancels F(Ei) F(E2). Our task consists in trans- 
forming this qualitative expectation into a proof. To this 
end we adopt the procedure of Ref. [] and simplify the 
term trg(QT 3 ) in Eq. (|J) in two steps. First, we write 
T as the product T = T C T D where T c (T D ) does not 
(does) commute with r 3 . With Qc — Tq LTc, this 
yields trg(Qr 3 ) — trg(QcT 3 ). Technically, the matrices 
Tc {Td) are obtained by exponentiating the coset gener- 
ators which anticommutc (commute) with t 3 . The result 
can be expressed in terms of the Ak x 4fc graded matrices 



Here, Q = T~ 1 LT is a graded matrix of dimen- 
sion 8k, d/i(T) is the invariant measure for inte- 
gration over the matrices T belonging to the coset 
U(2fc,2fc/4fc)/[U(2fc/2fc)] 2 , and trg denotes the graded 
trace. The 8k matrix indices are labelled prja which fol- 
low in lexicographical order. Here, p = 1, 2 distinguishes 
the retarded and advanced form of the Green functions, 
r = 1, 2 the energy arguments E± and E 2 , j — 1, . . . , k, 
and a is the supersymmetry index which assumes the 
valu es b for Bosons and / for Fermions. The matrices 
L and t 3 are diagonal matrices of dimension 8fc with en- 
tries ±1. In the matrix L (r 3 ), the two signs distinguish 
retarded and advanced Green functions (the two observ- 
ables F(Ei,H) and F(E 2 ,H), respectively). We have 
defined s = (E\ — E 2 )/2. Without loss of generality, we 
take £ > in the sequel. The explicit form of S depends 
on the form of the observable F(E, H) under considera- 
tion and is immaterial for what follows. Suffice it to say 
that S does not depend upon e, and that it contains a 
factor which depends exponentially on all the energy ar- 
guments LOj with j = 1 , . . . , k and a second factor which 
contains the source terms and, possibly, the coupling to 
open channels. 

Given Eq. (||) , the proof of ergodicity is a non-trivial 
and exact technical exercise which is sketched in the re- 
mainder of the paper. The limits of validity of the proof 
are those of Eq. (||) : We only admit ensembles of unitary 
matrices for which (i) the spectrum is confined to a sim- 
ply connected compact domain, and for which (ii) the 
fluctuation properties are translationally invariant (lo- 
cally stationary) in energy. Counterexamples to (i, ii) 
are soft confining potentials [] which typically have non- 
stationary fluctuation properties and for which another 
proof of ergodicity would be required. Counterexamples 
to (ii) are chiral random-matrix models where the point 
E = plays a distinct role. It is known that ergodicity 
does not hold near this point, (iii) Eq. (||) is restricted to 
the interior domain of the spectrum. Exclusion of the end 
points is due to a technical limitation of the method []. 
These points are of little interest for most applications 
of random-matrix theory where the energies E\ , E 2 are 
chosen to lie at or near the center of the spectrum. 

Eq. (||) looks like a Fourier integral. Given sufficient 
regularity of S near trg(Qr 3 ) = 0, we expect F(E 1 )F(E 2 ) 



T q = 



(4) 



with q = 1,...,4. These matrices belong to the coset 
U(k 1 k/2k)/[U(k/k)] 2 . The matrix elements of ij 2 and 
their conjugates in i 2 . 1 represent the "Cartesian" coor- 
dinates of T q . As usual, t q 2 and t 2 } are symmetry- 
related. Moreover, we have t^ 1 = (1 + t^t 11 ) 1 / 2 and 
i 22 = (1 + tft 12 ) 1 / 2 . The matrices T D (T c ) are given in 
terms of T\,T 2 (T 3 ,T4, respectively). Interest focuses on 
Tc which has the form 



A 11 



Tc = 





f n 

l 4 
'4 







z 12 

j-22 
'4 








j-22 
z 3 



(5) 



A second simplification of trg(QT 3 ) arises when, following 
Efetov [], we express T s for s = 3, 4 in "polar coordinates" 
by diagonalizing t 12 and t 21 . We write t 12 = UjA s (m 2 ) _1 
and t 21 = m 2 A s (uJ) _1 where the u p s with p = 1,2 are 
2kx2k graded matrices of the coset U(k/k)/[U(l)] 2k , and 
the X s and A s are diagonal 2k x 2k matrices related by 
symmetry. Then trg(Qr 3 ) = 4^ s ( — ) s trg(A s A s ) depends 
only upon the "eigenvalues" A s and A s of t 12 and t 21 , re- 
spectively. We write the eigenvalues in the form A^ = 

isin(0°./2)exp(i(%) and T sj = isin(6^/2) exp(-i<^.). 
The angle 9fj is positive imaginary (positive real) for 
a = b (a = /, respectively). The transformation is 
made unique by the requirement that for fixed s and 
a, the absolute values \8^\ decrease monotonically with 
increasing j. We write the matrices uf as products of 
two factors obtained by exponentiating the coset gen- 
erators diagonal (nondiagonal) in j. For the genera- 
tors diagonal in j, this yields the block-diagonal ma- 
trix iff = diag(w^), with w^j belonging to the coset 
U(1/1)/[U(1)] 2 and parametrized by two anticommut- 
ing variables 7^ , 7^* . By construction, the variables 
0"j , <j)"j , 7^ , and 7^* appear in Tc only in terms of the 
matrices \sj = w l SJ -K 3 {w^y 1 , where X sj = diag(A^), 
and in terms of the symmetry-related quantities Xsj ■ 

The Berezinian of the coordinate transformation to po- 
lar coordinates is singular. This is due to the occurrence 
of the factor ris=3,4 n.j'( cos ~ cos 

f 8j )~ 2 . Obviously, 
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with the assumed ordering of \0fj\, singularities occur 

when at least one of the two pairs 9 h sk , 6{. k , s = 3,4, van- 
ishes. This is what happens at trg(Qr3) = 0. The explicit 
form of the transformation i* 2 = li^A^u 2 ) -1 implies that 
the integral of S(Q) over the matrices u p s vanishes at the 
singularities, and is a regular function of the angles 6fj 
in their neighborhood. We handle the singularities in 
analogy to Refs. []. Lack of space forces us to sketch the 
essentials only, see Ref. []. 

We consider an integral X(kk) of the form of Eq. (||) 
but with an infinitesimal neighborhood of size 77 — » of 
9 h sk = 9^ k = 0, s = 3,4 removed from the domain of in- 
tegration. The removal is accomplished by inserting the 
Heaviside functions ]J S 9(\ sin(^ fc /2)| 2 + | sin(6»f fc /2) | 2 — r/) 
in the integrand. The resulting integral is evaluated in 
two ways. First, we show that the limit r/ — > exists. 
We do so by fixing the order of integrations in the in- 
tegral: We integrate over the angles Ofj after having 
integrated over all other variables. The result differs 
from the integral X = F(E 1 )F(E 2 ) in Eq. (§), how- 
ever. This is because when expressed in the original 
variables, the angles 0" fc possess nilpotent parts. Sec- 
ond, we evaluate the same integral by passing from the 
coordinates z s j — (9"j , 4>"j , 7^ , 7s 3 * ) to the coordinates 
x sj = (Xsj" i (Xs/* )*); where \1f are elements of the ma- 
trices \sj = w lj^sj(w 2 j)~ 1 introduced above. Since the 
Berezinian of this transformation removes the singularity 
from the integration measure, the limit can be taken by 
setting 77 equal to zero. Expanding the Heaviside func- 
tions in powers of the new anticommuting variables yields 
l(kk) as a sum of two terms. The first term is equal to 
the integral X. The second term contains delta functions 
of x s k in the integrand. These delta functions are gener- 
ated by the derivatives of the Heaviside functions. Thus, 
the second term is a "boundary" term where some of x s k 
are put equal to zero in the integrand and removed from 
the integration. We conclude that the integral X can 
be found by evaluating the contributions stemming from 
this boundary term and from the "volume" term X{kk). 
Evaluating the boundary term in polar coordinates, we 
encounter again a singularity. We handle this singularity 
in the same fashion as before. Continuing this proce- 
dure finally yields the integral F(E\)F(E2) as a sum of 
integrals over the polar coordinates generated in the de- 
scribed way. The sum contains the volume term and a 
number of boundary terms (Efetov-Wegner terms) of de- 
creasing complexity. The result for the autocorrelation 
function Cg reads 

fe / 

Cs = J2 T ^ ■ ( 6 ) 

h,U=0 

The prime on the summation sign indicates that the term 
with I3 = I4 = is omitted. This is because it is cancelled 
by the unlinked term F(E 1 ) F(E 2 ). The symbols X(l 3 l 4 ) 



stand for integrals of the form of Eq. (|3|), written in polar 
coordinates, with all variables z S j carrying the indices 
j > l a put equal to zero in the integrand, and removed 
from the integration. The integration over the angles 0^ 
must be done after all other integrals have been worked 
out. Thus, 

X(l 3 l 4 ) = J dfH ah {0)S hU (O) 

x exp(4^(e/d)^(-) s trgsin 2 (^/2)) , (7) 

sj 

where d/j; 3 ; 4 (#) denotes the measure for integration over 
the angles 0^-, and Si 3 i 4 (9) stands for the integral of S(Q) 
over all remaining variables. By construction, Sj 3 ; 4 is an 
even function of each 0"-, and, for fixed s and a, a sym- 
metric function of the 9% 's with j <l s . For 0^ , 0^ -> 
at fixed values of the remaining angles 0"-, S; 3 ; 4 tends to 

zero with the first power of \6 b sls | 2 + \0' sl | 2 . 

The entire e-dependence of the integrals X(l 3 l4) resides 
in the exponential in Eq. (Q). Therefore, the asymp- 
totic behavior ((e/d) — + 00) of these integrals is de- 
termined by the behavior of the integrands for small 
values of the 0™ 's. We introduce the rescaled angles 
6£j — (e/d) 1 / 2 9fj and expand all terms in the integrand 
in powers of (e/d) -1 / 2 . Since all terms are even functions 
of each 0"-, the expansions proceed in powers of (e/d)^ 1 . 

As usual, we extend the domain of integration over 0^ ■ 
to infinity. The asymptotic behavior of Z(/ 3 ?4) is given 
in terms of the leading terms in the expansions. For the 
measure and the exponential, these terms can be straight- 
forwardly found. The detailed form of the leading term 
in the expansion of 5 , ; 3 i 4 (0) is specific for the observable 
under consideration, and will change when another such 
function is considered. However, the general properties of 
Si 3 i 4 discussed above imply a lower bound on the degree 
ni 3 i i of this term. For Z3, 14 ^ 0, we find that ni 3 i 4 must 
be larger than or equal to 2(7 3 +Z 4 ) — 2; for l 3 — (I4 = 0), 
the corresponding bound is (2l 4 — 1) ( (2/ 3 — 1), respec- 
tively ). The asymptotic behavior of the autocorrelation 
function Cg is, therefore, determined by the terms X^^l^) 
with the lowest values of I3J4, (I3J4) — (1, 0)(0, 1)(1, 1). 
The function Cf is symmetric in Ei,E%, and only even 
powers of (e/d)' 1 can, therefore, occur in the asymptotic 
series. We conclude that the autocorrelation function 
Cs vanishes asymptotically at least as the second inverse 
power of (e/d), 

C £ - const x (e/d)' 2 . (8) 

This behavior is generic and guarantees ergodicity. 

We point out that nowhere in this proof have we used 
specific properties of the function S appearing under the 
integral in Eq. (H). Thus, our argument applies equally 
to all 4fc-point functions. But k is an arbitrary posi- 
tive integer. Hence, our result applies to all GUE ob- 
servables which depend on energy. We have, however, 
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used that F(E, H) is an observable, hence real, and that 
Cg is a symmetric function of E± , E% . For a non-real 
quantity G(E,H), the associated correlation function 
G £ = G{E 1 )G{E 2 )* - G(E^ G(E 2 )* is not symmetric 
in E\ , E% and generically decays asymptotically only like 
(e/cf) -1 . The S'-matrix autocorrelation function worked 
out in Ref. [] serves as an example. 

We turn to the second case and consider an observ- 
able F(B) = F(H(B)) which depends on the magnetic 
field strength B. The magnetoconductance autocorre- 
lation function is an example. The running average is 
performed over B. The arguments leading to Eq. (|^) 
apply analogously. Ergodicity now requires vanishing of 
the autocorrelation function of F(B) at large field differ- 
ences. A novel problem arises, however, in the modeling 
of the dependence of the random-matrix ensemble on B. 
For the case of the energy E, there is no such problem 
because the energy dependence of any observable arises 
from its dependence on a retarded or advanced Green 
function (_E ± — H)^ 1 . This leads to the e-dependence 
displayed in Eq. (||) and used throughout. In the case of 
a magnetic field, the Hamiltonian itself carries the para- 
metric £?-dependence. We again consider F(B) to be a 
fc-fold correlation function and ask: How do we model 
the random-matrix ensembles that describe the Hamil- 
tonians referring to 2k different magnetic field strengths 
B ± AB + f3j, j = 1, . . . , fc? Here, B is some mean mag- 
netic field strength, the variables B ± AB correspond 
to the variables E\ and E 2 in the case of an energy- 
dependent observable, and the variables f3j are analogous 
to the uij's. We consider the case where the autocorre- 
lation function is invariant with regard to reversing the 
sign of each of the (3j's, and model the 2k Hamiltonians 
by the matrices 

HftP = H$ + (-Y+XVWH^ + b jy /T/NHW . 

(9) 

Here, bo is proportional to AB, and the bj's with j = 
1, ... , k are proportional to the fij's, respectively. The 
index r takes the values r = 1, 2 for B + AB and B — 
AB, respectively. The three GUE ensembles H"' with 
I = I, II, III are taken to be uncorrelated. Eq. (0) is a 
straightforward generalization of the usual model applied 
in the case of two different magnetic field strengths, see 
Ref. 0- 

We consider the limit \AB\ — > oo of the au- 
t ocorrelation function C B = F(B + AB)F(B - AB) 
—F(B + AB) F(B - AB), keeping the 0/s fixed. To 
calculate Cb, we again use the supersymmetry technique 
and arrive at an expression which has the form of Eq. (|^) 
except that the exponential function under the integral is 
replaced by exp(— [6q/2] trg(Qr3) 2 ). Moreover, the func- 
tion S now depends on the field variables bj,j = 1, . .., k 
rather than on the energies u>j . We proceed as before and 
arrive at an equation which has the form of Eq. (p) except 



that now the integrals 1(13,14) contain the exponential 
exp(26p trg sin 2 6 S j). This expression is a symmet- 
ric function of AB and, thus, automatically obeys the 
symmetry property of Cb- Proceeding as before we fi- 
nally conclude that the leading term in the asymptotic 
series for Cb stems generically from the Efetov-Wegner 
terms X(10) and 1(01) and is of order b^ 2 . This implies 
that for I bo | — ► 00, we get 

Cb ~ const x b^ 2 . (10) 

This asymptotic behavior is sufficient to guarantee er- 
godicity. 

In summary, we have shown ergodicity to hold for a 
wide class of ensembles of unitary random matrices and 
for all observables which depend either on energy or mag- 
netic field strength. We expect that our proof applies 
likewise to the orthogonal and to the symplectic ensem- 
bles. 
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in Heidelberg for their hospitality and support, and ac- 
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GAUK in Prague. 
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